Critical conductance of a one-dimensional doped Mott insulator 
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We consider the two-terminal conductance of a one-dimensional Mott insulator undergoing the 
commensurate-incommensurate quantum phase transition to a conducting state. We treat the leads 
as Luttinger liquids. At a specific value of compressibility of the leads, corresponding to the Luther- 
Emery point, the conductance can be described in terms of the free propagation of non-interacting 
fermions with charge e/v2. At that point, the temperature dependence of the conductance across 
the quantum phase transition is described by a Fermi function. The deviation from the Luther- 
Emery point in the leads changes the temperature dependence qualitatively. In the metallic state, 
the low-temperature conductance is determined by the properties of the leads, and is described by 
the conventional Luttinger-liquid theory. In the insulating state, conductance occurs via activation 
of e/\/2 charges, and is independent of the Luttinger-liquid compressibility. 

PACS numbers: 71.27.+a 71.10.-w 73.43.Nq 73.21.Hb 
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I. INTRODUCTION AND OUTLINE 

The dc conductance of a one-dimensional (ID) system 
crucially depends on the properties of the leads attached 
to it. Unlike in higher dimensions, this dependence is 
not merely a technological, but rather a generic physics 
problem. 

In the case of the conventional Luttinger-liquid model 
for a quantum wireji the role of the leads in the zero- 
temperature (T = 0) conductance is well studiedii^^^ 
As long as the electrons do not backscatter off inhomo- 
geneities, the conductance is not sensitive to the nature 
of carriers in the wire, and is determined solely by the 
leads. In the simplest case of the noninteracting leads, 
the conductance is quantized in units of e 2 /27r?i per chan- 
nel (each of the two spin polarizations counts for a sepa- 
rate channel). This conclusion also remains true for the 
incommensurate phase of a so-called Mott wire (doped 
ID Mott insulator)^ 7 . Thus the charge fractionalization 
characteristic of these systems^ is not observable in the 
zero-temperature ballistic conductance; if the leads are 
Fermi liquids, then the ballistic conductance of the wire 
appears to coincide with that of the free electrons^ 

The Tomonaga-Luttinger model is applicable to a rel- 
atively dense system of particles forming a compressible 
liquid. The "distance" from the incompressible state may 
be characterized by a parameter — r, having the meaning 
of chemical potential measured from the point where the 
system becomes incompressible. For example, for free 
fermions filling up a conduction band 



— r = /i — A 



(1.1) 



with [i and A being the Fermi energy and band edge, 
respectively. At T = and [i < A the system is incom- 
pressible. 

The conventional Luttinger liquid model requires rela- 
tively high carrier density, — r/T 1. How robust is the 



apparent free-electron picture of transport if the condi- 
tion —r/T 3> 1 is violated? One may expect deviations 
from the apparent free-fermion form of the two-terminal 
conductance in some intermediate region of parameters?: 
in a conventional quantum wire. However, for fully spin- 
polarized electrons (equivalent to spinless fermions) at 
the point r — of the quantum phase transition the in- 
teraction becomes irrelevant. The corresponding conduc- 
tance equals | x e 2 /h, i.e. a half of the unit conductance 
quantum e 2 /h for spinless fermions, which is again indis- 
tingushable from a free-fermion result ~ 

Charge carriers at the commensurate-incommensurate 
quantum phase transition in a Mott wire can also be 
mapped onto free spinless fermions ~ The distance to the 
phase transition is characterized by a chemical poten- 
tial — r, see Eq. (jl.ip . here measured with respect to the 
Mott gap, A. However, these effective fermions have frac- 
tional charge, and are related to the original electrons in a 
complicated way. It is therefore an important unresolved 
question as to what the two-terminal critical conductance 
of a Mott wire is. 

The aim of this work is to investigate the finite- 
temperature conductance around the quantum phase 
transition point in a Mott wire. Our main finding is 
that such a setup may provide the way to experimentally 
access the elusive character of the fractional quasiparticle 
charge in an interacting ID system. 

To be specific, we consider the two-terminal conduc- 
tance of a one-dimensional electron system close to the 
Mott transition. The details of how the underlying com- 
mensuration is realized in practice, e.g. by a periodic po- 
tential, will not be important; the only parameter from 
the Mott wire we will need in the end is the effective value 
A of the charge gap. Next, we assume that this wire is 
smoothly attached to uniform identical Luttinger liquid 
leads characterized by the Luttinger parameter Kr,. We 
evaluate the finite-temperature conductance close to the 
Mott transition perturbatively in Kr, — i. The value of 
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Kl = h is special (Luther-Emery poinfei£), as it allows a 
mapping of the entire system onto a free-fermion one. In 
lowest order in Kl 



|, we find 



A> 



r - 2e 



(1.2) 



where r is the tuning parameter (jl.l[) of the Mott tran- 
sition, Gq is the conductance quantum, and 



finding its two-terminal dc conductance. The latter is 
affected by the different nature of the charge carriers in 
the leads and in the wire. 

The paper is organized as follows. In Section [IT] we 
consider the properties of the uniform Mott wire. In Sec- 
tion IIIII we consider the Mott wire with the leads, and 
evaluate the modification of the critical conductance in- 
duced by the presence of the leads. Section HVl is devoted 
to the discussion of our results. 



(1.3) 



is the Fermi function. On the metallic side, r < 0, the 
conductance approaches the value KlGq in the limit of 
zero temperature, i.e., it is determined by the Luttinger 
liquid parameter of the leads in agreement with Ref. |g. 
On the other hand, on the insulating side of the transi- 
tion, r > 0, transport at low T is thermally activated 
as the Fermi function reduces to a Boltzmann factor, 
f(r/T) e~ r l T . In lowest order in the fugacity e~ r l T , 
we find that the conductance G — ^Goe~ r ^ T is not af- 
fected by the leads, i.e., is independent of Kl- In par- 
ticular, the prefactor of the exponential is determined by 
the effective conductance quantum \G$ attributed to the 
degrees of freedom of the critical Mott wire. As we will 
explain in detail below, the critical degrees of freedom 
are the fermions that carry fractional charge 



eci 



(1.4) 



At low temperatures transport in the Mott insulator oc- 
curs via activation of the fractional charges (|1.4p , result- 
ing in the announced prefactor \Gq = 2e^j l /h. Finally, 
in the extended quantum critical regime, \r\ <C T, the 
conductance is dominated by its critical value, 



Gcr = G\ r _ — 



2K L + 1 



G + O{K L ~\f 



(1.5) 



Note that the specific fraction of the conductance quan- 
tum here depends both on the fractional charge of the 
carriers and on the properties of the leads. These results 
differ from the considerations of Mori et al.£ 

Our study of the two-terminal conductance of the criti- 
cal Mott insulator complements the theoretical investiga- 
tion of ballistic transport in massive theories in general, 
for a recent review see Ref. [l3|. The thermally activated 
quasiparticles above a gapped ground state in a homoge- 
neous system (no leads) may travel ballistically thus giv- 
ing rise to a singular contribution to the optical conduc- 
tivity, Using, (w) = 2icD8(u), with a characteristic Drude 
weight D = D(T). The low-temperature Drude weight of 
the Mott insulator is believed to remain finite even away 
from criticality^ (In the exact commensurate case of an 
undoped Mott insulator the fate of the Drude weight is 
less clear and remains a subject of ongoing research i 15 i 16 ) 
Remarkably, the knowledge of the Drude weight in the 
conductivity of a uniform ID system is insufficient for 



II. UNIFORM MOTT WIRE 

In this Section we consider the model of the uniform ID 
Mott wire. In Sec. Ill Al we introduce the Hamiltonian of 
our system, and in Sec. Ill Bl we sketch the critical theory 
describing the Mott transition. In Sec. Ill CI we outline 
the application of the Kubo formula for the calculation 
of conductance, and derive the critical conductance for 
the Mott wire without leads in Sec. Ill Dl 



A. Model of the one-dimensional Mott insulator 

Consider an interacting one-dimensional electron sys- 
tem in the vicinity of the Mott transition. The latter can 
be realized, e.g., in a Hubbard model close to half- filling. 
The effective Hamiltonian is conventionally represented 
in terms of the bosonic charge and spin excitations that 
decouple at sufficiently low energies £ In this work we 
only focus on the charge sector which, near the transi- 
tion, is described by the sine-Gordon model 8 



n[<t> c , e c 



dx 
2^ 



Ku(d x e c ) 2 + -{d x <p c f 



(2.1) 



Here K is the Luttinger parameter, u is the plasmon 
velocity, Vumki is the Umklapp scattering amplitude 
(given by the Fourier component of the lattice poten- 
tial at the reciprocal lattice vector) that corresponds 
to (? 3 in the standard <?-ology classification^ and a is 
a short distance cutoff. The conjugate fields <p c and 
6 C , with \~d x (j) c {x), c {x')] = —iS(x — x'), describe the 
smooth components of the density and current fluctua- 
tions within the system, 



Sn = -V2 



d x c. 



and j = V2K: 



(2.2) 



respectively. The prefactor v2 in these expressions is at- 
tributed to the two spin polarizations of the electrons. 
The parameter go measures the distance to half-filling. 
In the presence of repulsive interaction between the elec- 
tons, K < 1, the umklapp scattering process, Vumki, be- 
comes relevant and opens up a gap in the spectrum at 
half filling qo = 0, rendering the system into the Mott 
insulating phase. 
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The elementary excitations of the Mott insulator are 
massive quantum solitons (kinks) of the model (|2.ip . 
Classically, the kink must connect two discrete values of 
4> c {x — > ±00), which realize the minimum of the cosine 
term. This dictates the fractional value of the topological 
charge 



Q = i 

7T 



dx d T 



olito 



V2 



(2.3) 



Below we will absorb the topological charge (|2.3|) carried 
by the solitons into the fractional electric charge quan- 
tum, e C i = Qe, cf. Eq. (fi~4|) . 

Upon detuning the system from half-filling by in- 
creasing qo across a critical value the system becomes 
conducting via the commensurate-incommensurate (CI) 
transitions^ The ground state of this conducting phase 
is characterized by a finite topological charge density. 
Consequently, this transition is described 18 ' 19 by exploit- 
ing the duality 21 ^! between the sine-Gordon model (|2.1j) 
and the massive Thirring model. To establish this dual- 
ity, we make a canonical transformation to the rescaled 
fields 

<j> = V20 c , and 6 = 6 C /V2. (2.4) 
In these variables the model (12.111 becomes 



dx r 

27 



(2.5) 



2Ku{d x 6f + — (d x( f) 2 

, 2Vumkl ,„ , \ 

dx 7K — W cos(20 - q x) . 



The above Hamiltonian for the charge sector can now be 
refermionized with the help of the standard bosonization 
formula. 



1 



(2.6) 



with k = +1 and —1 for right (R) and left (L) movers, 
respectively. The fermionic particles associated with the 
operator ip K are identified with the solitonic excitations 
of the sine-Gordon theory; the fermionic density thereby 
coincides with their topological charge density. The re- 
sulting fermionic theory is the massive Thirring model, 

Hm = J dx {V [-iva 3 d x -fx + a 1 A] * (2.7) 

+| (.94 + 92) (*^) 2 + \ (54 - 92) (*V 3 *) 2 } 

where a 1 are the Pauli matrices, and we used the spinor 
notation \ft = (^p, ip L ) for the new spinless fermions. As 
the Fermi momentum for the fictitious Dirac fermions \& 
is set to zero (at the Dirac point), there are no spatially 
oscillating factors in Eq. (|2.6[) . These fermions (or the 
solitonic excitations of the sine-Gordon theory) describe 
the density excitations smooth on the scale of the original 
lattice responsible for the Mott phase. The parameters 



of the model (|2.7p are the chemical potential /i = —vqo/2 
and the gap A = Vijmki/(27ra)j22 the velocity v and the 
interaction constants g± and 91 are implicitly given by 



u = v\J (1 + vg 4 f - {vg 2 y 



2K 



1 + vgi - vg 2 



1 + vgA + vgi ' 
where the density of states for the ^/-fermions 

1 

2nv 



(2.8) 
(2.9) 

(2.10) 



Note the redundancy of the g^ interaction as it enters the 
mapping only in combination with the velocity, 2irv-\-g±. 
We deliberately choose to keep both the 52 and g± in- 
teraction constants, as it will help us later to classify 
different perturbative contributions to the conductivity. 
The Luttinger interaction parameter for the new fermions 
(|2.7p is 2K. As a result, the interaction among them 
vanishes at the Luther-Emery point, 12 K = i. Later, 
in Section IIIII we will apply the above formalism to the 
wire with the leads, and develop a perturbative expan- 
sion around this special value for the Luttinger parameter 
within the leads, Kl. 

In terms of the spinless fermionic degrees of freedom 
the charge current density, J = ej, is given by 



d 9 

J = \T2eKu-^- = V2e C i^vja 3 ^. 



(2.11) 



Again, the overall multiplicative factor v2 can be traced 
back to the spin degree of freedom of the original elec- 
trons. As announced, we absorbed the topological charge 
associated with the fermions, Q = 1/V2, into the frac- 
tional electric charge quantum, eci = Qe, of Eq. (|1.4[) . 
Moreover, the charge velocity of the new fermions reads 



vj = v (1 + vgi - vgi) 



(2.12) 



For the general case 32 7^J74> the velocity vj depends on 
the interactions (cf. Ref. [g, Sec. 7.2). 

The quadratic part of the refermionized theory (|2.7|) 
can be diagonalized by a Bogoliubov transformation re- 
sulting in the two bands with the spectrum ±e^, where 
€k = -\/A 2 + (vk) 2 . When the chemical potential lies 
within the gap, < A, the ground state is a Mott in- 
sulator. When the chemical potential is tuned into one 
of the two bands the ground state becomes a conductor 
via the CI transition. 



B. Critical theory of the 
commensurate-incommensurate Mott transition 

From now on in this work we will focus on the prop- 
erties of the system close to the Mott transition and as- 
sume in the following that the chemical potential of the 
\1/ fermions is near the upper band edge, fj, « A, i.e., 
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r « 0, where r is the tuning parameter for the quantum 
phase transition, cf. Eq. The generalizaton to the 

situation when fi w —A is straightforward and will not 
be further discussed. 

The effective Hamiltonian governing the transition is 
given by nonrelativistic free spinless fcrmionS ) 18 ' 19 



dx (x) 



2m 



c(x), 



(2.13) 



where the operator c 1 creates a fermionic excitation in the 
upper band, the mass is m — A/v 2 , and r is defined in 
Eq. (jl.ip . The scaling dimension of the tuning parameter 
and the critical dynamical exponent are dim [r] = \/v r = 
2 and z = 2, respectively. 

As the gas of fermions is very dilute at the transition, 
their interaction is weak. As the CI transition is ap- 
proached, r — > 0, the interaction among the fermions de- 
creases faster than their density, rendering the c fermions 
free at the transition. Speaking crudely, the fermions are 
free not because they do not interact per se but because 
there are too few particles to interact with. Formally, this 
can be seen from considering the residual interaction 



n/CI 
ri 'mt 



= 7 6k + k '>P 



P+P' ^kk';pp' C k C k' C p C p' 



t J 



kk' pp' 



The interaction amplitude Tkk' ;pp' is strongly constrained 
as the Pauli principle requires the effective amplitude to 
be antisymmetric upon exchange of momenta p <-> p' and 
k <-> k'. Near the CI transition we can expand the am- 
plitude, and the lowest order contribution is suppressed 
by two powers of momenta, 



kk' \pp' 



= 1 (k-k')(p-p'), 



(2.15) 



with the coefficient 7 = — g2V 2 / (4A 2 ). As a consequence, 
at the quantum critical point the residual interaction, 7, 
is irrelevant in the renormalization group sense. 

The phase diagram of the CI Mott transition is shown 
in Fig. [T] The tuning parameter r controls the distance 
to the quantum critical point; for r > the system is in 
the commensurate phase with an insulating ground state. 
In the metallic phase, r < 0, at lowest temperatures 
T < III, with log(|r|/T L L) oc l/|r|, a crossover between 
the critical behavior and the Luttinger liquid takes place; 
the residual interaction 7 at finite density induces Lut- 
tinger liquid correlations in the spectral function of the c 
fermions with the Luttinger parameter K c — 1 cx 



C. Transport 

We will consider electric transport through the criti- 
cal Mott wire. Like the thermodynamic properties, the 
transport has universal characteristics at the Mott tran- 
sition, on which we focus in the following. The conduc- 
tance can be obtained from the zero frequency limit of 



T . 
A - 



quantum critical ' 

regime /t ~ r UrZ 



»T LL 



\ t 







FIG. 1: Phase diagram of the CI transition for temperatures 
smaller than the gap A, see text. A quantum phase transition 
with exponents v r = 1/2 and z — 2 separates the commen- 
surate, r > 0, and the incommensurate, r < 0, phases. The 
CI Hamiltonian (12.13[) controls the vicinity of the quantum 
critical point. The quantum critical regime occurs within a 
coneT~ \r\ v ^ = \r\. 



(2.14) tue conductivity 



23 



G = lim a(x, x ; u>), 



(2.16) 



at fixed positions x and x' . The conductivity a(x, x'; uS) 
describes the linear response of the charge current at po- 
sition x upon applying an external— electric field £(x' , u>) 
at position x' with frequency oj 



(J(x,cu)) = / dx'cr(x,x';uj)£(x' ,uj). 



(2.17) 



The conductivity is derived with the help of the Kubo 
formula that relates a to the charge current autocorrela- 
tion function 



a(x, x';lu) — 



(2.18) 



where the retarded function is defined by 



' / dte^{[J{x,tl J(x',0)}). (2.19) 



The subtraction of the uj = mode in (|2.18p ensures 
gauge invariance. In the following, it will be convenient 
to obtain the retarded response function from the tem- 
perature correlation function 

ICx,x'-n n =f f dTe innT (T r J(x,T)J(x',0)) (2.20) 
Jo 

by the analytic continuation from the upper half-plane, 
fc r x%'-,u = £x,x';Si n \in n ->u;+io, where Q n = 27mT is the 
bosonic Matsubara frequency. 
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D. Critical conductance of the homogeneous 
system 

Consider first the uniform system close to the metal- 
insulator transition. At the criticality, the transport is 
determined by the free nonrelativistic c-fermions, i.e., by 
the Hamiltonian (|2 . 13[) . The resulting conductance 



properties of the critical conductance are modified. Be- 
low we develop the perturbation theory in the deviation 
from the Luther-Emery point, and evaluate the critical 
conductance in the first order in Kl — i. 



G (h 



<HU) = ^f(r/T) 



1 



G f(r/T) 



(2.21) 



A. Leads with Luttinger parameter Kl 



is simply the conductance quantum 2eQ l /h for the ef- 
fective fermions with charge (jl.4| . times the occupation 
factor that is given by the Fermi distribution (|1.3[) . 

The critical conductance is universal and its tempera- 
ture dependence is controlled by a Fermi function (|1.3p . 
The non-universal corrections arising from the residual 
interaction (|2.14p vanish with the distance to the quan- 
tum critical point as the interaction 7 is an irrelevant 
perturbation. In particular, in the limit of low temper- 
atures the conductance at criticality, r — 0, approaches 
the value 



^(ho 



— -7G0 

r=0,T^Q 4 



(2.22) 



T=0,r^0+ 



At zero temperature the conductance G^ ho 
in the insulating phase, while in the metallic phase 
=o,r^o_ = 9 Go, thus the conductance exhibits 
a jump across the transition of the height 



G (hc 



T=0.r^0- 



1 



-G . (2.23) 



The factor ^ originates from the fractional nature of the 
charge carried by the critical fermionic degrees of free- 
dom, eci = e/ \/2. As was pointed out in Ref . [25|, the 
conductance jump across the transition (|2.23[) in fact cor- 
responds to a full conductance quantum but with frac- 
tional charge, Go/2 = 2e 2 cl jh. 

In the following section the modifications to the critical 
conductance (|2.21[) due to the presence of leads attached 
to the Mott wire are discussed. 



The fermions entering the Thirring model (|2.7p are 
the appropriate degrees of freedom describing the Mott 
transition as the interaction, gi and (74, among them is 
irrelevant. In the following we will also adopt a descrip- 
tion of the leads, where the Mott gap is absent, A = 0, 
in terms of the same \l/-fermions. The advantage of such 
a description is that the Mott region of the wire can be 
incorporated into an effective scattering approach for the 

fermions which is outlined in the following. The caveat 
is that, unfortunately, these effective fermions will in gen- 
eral, i.e. for Kl ^ i, interact within the leads. This is 
what makes this problem difficult, as one has to solve the 
strongly interacting fermionic problem which also lacks 
translation invariance. The influence of the interaction 
within the leads is considered in Section IIII Bl perturba- 
tively. 

Mori et al. in Ref. iTj considered a very similar setup 
as we do here. They also calculated the current-current 
correlation function with the help of the \P fermions at 
the Luther- Emery point but only for the uniform system, 
where the Mott gap extends over the full wire. They 
proceeded by using the expression for the homogeneous 
case to derive the conductance in the presence of the leads 
where the original electrons do not interact. We think 
that the last step has led them to the erroneous result for 
the temperature dependence of the critical conductance. 
We will further elaborate on the difference between the 
treatment of Mori et al. and ours in Section IIV1 



III. CRITICAL CONDUCTANCE IN THE 
PRESENCE OF LEADS 

In a two-terminal conductance measurement the crit- 
ical Mott wire will be connected to leads. In the pres- 
ence of leads the conductance will differ from the ex- 
pression (|2.2ip for the homogeneous system. We assume 
here that the leads contain Luttinger liquids character- 
ized by the Luttinger parameter Kl ■ Here we argue that 
the critical conductance is that of the homogeneous case, 
Eq. (I2.2ip . only if the Luttinger parameter Kl has the 
Luther-Emery value Kl = 1/2. Note that this value of 
Kl corresponds to strong interaction between the original 
electrons in the leads. If Kl differs from \ , the universal 



1. Scattering states 

Let us consider the situation where the Luttinger pa- 
rameter of the leads Kl = \, corresponding to g 2 = 
in the Hamiltonian (|2.7p . Here we also set 34 = 0. The 
absence of interactions allows us to describe the wire in 
terms of a simple ID scattering problem. The leads can 
then be characterized by the states of the 'J fermions that 
scatter off the Mott region. The crucial point is that the 
latter will be treated as an effective point scatterer. This 
is a correct approximation for the dc Landauer trans- 
port, and physically it means that we attach infinitely 
long leads to the Mott region. 

The scattering wave incident from the left has the stan- 
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dard form 
1 



(3.1a) 



iex jv 



iex jv 



t L (e) 



lx l v r L u ) ( j) , x < 
x > 0. 



The energy e is measured here (and everywhere below) 
with respect to the Fermi level. The phase velocity in the 
R- and L- moving states (upper and lower spinor com- 
ponents) equals ±v correspondingly. The wave incident 
from the right reads 



—iexjv 



««w (?) 



(3.1b) 



'2,7TV | g—iex/v 



c/v 



r R {e) 



x < 
x > 0. 



The transmission and reflection coefficients are compo- 
nents of the scattering S-matrix, 



S(e) = 



_ ft L (e) r R (e) 



r L (e) t R (e) 



(3.2) 



which is unitary, S+S = 1. The unitarity of the S-matrix 
ensures that it has only three independent parameters, 
which can be chosen in the form 



ti 



tj 



rL 



J<P- 



re 



(3.3) 



with <t>t — {<t>+ + (f>-)/2, and positive r and t, r 2 + 1 2 = 1. 
The relation = t R is a consequence of the time-reversal 
symmetry. 

The presence of the Mott gap at the center of the 
wire makes the transmission and the reflection coeffi- 
cients strongly energy-dependent: The propagation is 
blocked for the fermionic states with energies within the 
gap, e < r. These states therefore have zero transmis- 
sion coefficient as we assume the Mott gapped region 
of the wire to be sufficiently long such that the tunnel- 
ing across it can be neglected. On the other hand, the 
fermionic states with energies above the gap, e > r, are 
able to traverse the Mott region. We will assume that 
the leads are attached to the Mott-gapped region of the 
wire in the adiabatic fashion (i.e. Vumki = Vumki(aO is a 
smooth function), such that the backscattering of parti- 
cles moving above the gap can be neglected. The tuning 
parameter, r, of the transition thus divides the fermionic 
scattering states into the propagating waves with unit 
transmission for energies e > r, and into the standing 
waves for energies e < r. As a result, we model the Mott 
region of the wire by utilizing simple unit step-function 
transmission and reflection amplitudes 



t( e ) = 0(e-r), r(e) = 6(r - e) . 
Summarizing, 

(C) ~ I -0(r - e)e^-^ 6(e - r)e*+*& 



(3.4) 



(3.5) 



The scattering phases 4>±(e) are accumulated while 
traversing or reflecting from the Mott barrier. The bar- 
rier shape is encoded in their smooth energy dependen- 
cies, 



4>±{t)~4P ± +l±e/v, 



(3.6) 



The scattering lengths £± are a measure of the effective 
length of the region of the wire where the Mott gap is 
present. 



2. Green's Function 

The single-particle Green's function in the Matsub- 
ara representation is a 2 x 2 matrix in the spinor space 
[a, 6=1,2 are the spinor components of the wave func- 
tions f3~T]) ] 

Gtx';r = - (Tr^ix, r)^ b (x', 0)) . (3.7a) 

In what follows we will be using its Lehmann represen- 
tation 



Gx,x' \u 



de- 



(3.7b) 



where A&,. e = £ ^A x )^ x ') > ( 3 - 7c ) 
<j=± 

and uj n = n(2n+ 1)T is a fermionic Matsubara frequency. 

As we show now, the Green's function (|3.7|) has quali- 
tatively different properties when its arguments x and x' 
are located on the same and on the opposite sides of the 
impermeable potential barrier. Introducing the labels 

s= signx and s = signs', (3-8) 

one finds: 
s + s' : 



D ie(x — x')/v-\-is<j>t (e) 

o 



o 

~ie(x — x )/v—is<fit( i e) 



(3.9a) 



gie(x— x ) I v 

o 





-ie(x — x )/v 



+ VS Q r -e 







ie(x+x')/v+is<fi s (e) 



-ie(x-\-x )/v—is4> s (e) 



o 



(3.9b) 



The spectral function (|3.9ap . as well as the first term in 
the spectral function (|3.9b[) become translation-invariant 
in the limit \x — x'\ — > oo, since they both depend only 
on the coordinate difference x — x' . They originate from 
the propagating states and from the corresponding co- 
moving parts of the standing waves. The second term 
in Eq. (|3.9b[) is due to the counter-moving parts of the 
standing waves and comes only from the states with en- 
ergies below the barrier, e < r. It lacks translational 
invariance, and is affected by the barrier shape through 
the energy-dependent reflection phase shifts <p±. 
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3. Polarization operator 

A basic ingredient in the following calculation is the 
polarization operator that is a convolution of the tem- 
perature Green functions (|3.7p . 

n';.', M; = -tJ2^- {r l g(x,x';uj n + n n yg(x' ,x-,uj n )} 

(3.10) 



where i, j = 0, 1 and 
'1 



T = U p 



with /3 l 



1 . i = 

-1, i: 



l| (3.H) 



The polarization matrix 11^ ,.^ depends on the bosonic 
Matsubara frequency f2„ = 2imT. A detailed calculation 
of the above expression can be found in the Appendix [A] 
In the dc limit (iCl n — > w + iO, lo — > 0) it reduces to 



n^ ; n„=n^„ + [/(r/T)-i] 



(3.12) 



1 - ss' 



-n 



(0) ij 



where / is the Fermi function (|1.3p . and we used the 
abbreviation (|3.8p . The polarization operator matrix in 
the absence of the Mott barrier is 



(0) ij 



O n ' 



= -27TZ/V 



| O n | |x — x' | /tJ 



(3.13) 
(3.14) 



IO, 



Vt n sign (x — x 1 ) 



Cl n sign (x — x') 



vanishes as sm(ux'/v) in the dc limit: Given enough 
time, the system finds out about the barrier. For finite 
/, the part of the expression (|3.12[) responsible for the 
transmission of the particles is proportional to f(r/T). 

Technically, the presence of the barrier leads to an ad- 
ditional term in the expression (|3.12[) that is multiplied 
by the factor f(r/T) — 1. This term can be further sim- 
plified by noting that, in the dc limit, the point-scatterer 
approximation for the Mott region allows us to substitute 
|x + x'| = |x| + |x'| for s = s', and |x — x'\ = \x\ + \x'\ for 
s s' , yielding 



irl°^ n = niV n e -l^'l/-, s ^s', (3.15) 



n 



(0)« 

X, — X ! \Vl-n 



- n tj 



\Q n x \/v 



(3.16) 



Using the above simplification the expression Eq. (|3.12[) 
can be cast into a form characteristic for a scattering 
problem, 



n, 



n 



(0) 



+ IT 



,0;fi< 



IE 



0.x';O„ 



(3.17) 



In the limit of small Q n (in the sense of analytic con- 
tinuation itt n — > lu — > specified above), the scattering 
T-matrix is given by 



1 = 




1 



T = 



1 - f(r/T) 

27r^ 2 |f7 T J 



(3.18) 



4- Critical conductance with leads at the Luther-Emery 
point Kl = | 



with the density of states (|2.10|) . The rows and columns 
of the matrix (|3.14p are numbered by the indices i,j = 
0, 1, such that fl 00 corresponds to the upper left cor- 
ner. The corrections to the result (|3.12p are subleading 
in powers of frequency fl n , see Appendix [Al for details. 

Let us discuss the meaning of the expression (|3.12D . 
The trivial observation is that, in the absence of the 
barrier, / = 1, and II = II^ ^ becomes translation in- 
variant. In the opposite case of an infinitely high bar- 
rier, / = 0, and n^/.jj = identically vanishes when 
the points x and x' are on the opposite sides of the 
barrier (no transmission). We now consider the polar- 
ization operator when x and x' are on the same side, 



n 



n 



(0)« 
x,x'\Q n 



This expression has 



a static part (that remains finite in the limit f2„ — » 0), 
and the dynamic part cx VL n . The static part of the com- 
pressibility n° x ,. On=0 = 2v8{x — x') is not affected by 
the barrier at all (this in fact is true for any value of 
/) , as the density of states in the leads is independent of 
the barrier. The dynamic part is more interesting. For 
/ = 0, it is given by ^ + i-W^ . This ex- 
pression can be understood by looking at the dynamic 
part of the current-current correlator, corresponding to 
i = j = 1. After the analytic continuation it behaves 



oc toe 



iux j v 



sin(wx'/u). Therefore, the conductivity (|2.18p 



From the evaluated polarization matrix (]3 . 1 7[) we can 
derive the conductance in the presence of leads that have 
a Luttinger parameter of the Luther-Emery value Kl = 
i The current autocorrelator (|2.20p is then given by 



r LE 



n 2 2tt11 



(3.19) 



The conductance f|2 . 16[) follows from the Kubo formula 



Gle- -G Q f(r/T), 



(3.20) 



with the conductance quantum Go = 2e 2 /h. We thus 
reproduce the critical conductance of the homogeneous 
system (|2.2ip . Below, we show how the critical proper- 
ties of the conductance are modified in first order in the 
deviation Kl — h- 



B. Critical conductance: Perturbation theory 
around the Luther-Emery point 

Here we evaluate the correction to the conductance 
arising from the deviation of the Luttinger parameter in 
the leads from its Luther-Emery value Kl = \- The 
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(a) 




(e) 



FIG. 2: First order corrections to the current autocorrelator. 
The straight line represent fermionic ^ propagators, the wig- 
gled line is the interaction and the dot is the current vertex 
^2e C i 



deviation from Ki, = | implies according to Eq. (|2.9[) a 
finite 172 interaction between the scattering \P fermions. 
In the following we calculate the conductance up to the 
first order in the interactions g 2 and 174. As in the first 
order the c/4 interaction only modifies the velocity, it is 
expected that it does not enter the conductance correc- 
tion. The verification of this expectation serves as a check 
of our calculation and is actually the only reason to intro- 
duce at all the redundant (74 interaction in the mapping 
(E^D and (EU). 

There are various contributions to conductance arising 
from finite g 2 and 34. First of all, these interactions mod- 
ify the expression of the charge current velocity vj, see 
Eq. (I2.12|) . As a consequence, the conductance obtains 
the modification, 



SG A /G LE = 2(g^~g 2 iy). 



(3.21) 



In addition, there are in total six first order diagrams 
that contribute to the current autocorrelator, see Fig. [5] 
Diagrams (a) and (b) cancel each other exactly due to the 
Fermi statistics. Diagram (c) does not modify the con- 
ductance because its contribution to o{x,x']uj) vanishes 
in the low frequency limit. In the absence of the barrier 
this is self-evident, as the small frequencies correspond to 
small incoming momenta that cannot change an R-mover 
into a L-mover. In the presence of the barrier this argu- 
ment essentially goes through, since only the propagat- 
ing parts (|3.9a[) of the Green's functions contribute; the 
latter have the same diagonal form as in the translation- 
invariant problem. Diagram (d) corresponds to the con- 
ductance correction coming from fermion scattering off 
a Friedel oscillation| 26 i 27 i 28 we will see that in our case 
this contribution is absent. The only modifications to the 
conductance arise from the self-energy diagram (e) and 
from the random-phase approximation (RPA) diagram 
(f). They are discussed in detail below. 



1. Scattering off Friedel oscillations: Diagram, (d) 

The self-energy correction in diagram (d) embodies the 
scattering off a Friedel oscillation. In the general case of 
a semitransparent barrier, this process leads to a loga- 
rithmic renormalization of the scattering coefficients and 
of the conductance ! 26 ! 27 ' 28 However, the S-matrix in our 
case has a very specific energy dependence (|3.5p as the 
transmission at a given energy is either zero or unity. In 
such a situation the renormalization e.g. of the trans- 
mission coefficient vanishes as it is proportional to the 
product of transmission and reflection coefficients at a 
given energy. As a consequence, diagram (d) does not 
contribute to the conductance here. 

The fact that there is no contribution from the Friedel 
oscillation in the first order in 172 can be undestood al- 
ready from the following simple argument. Consider the 
propagating state, for which the transmission is unity in 
the absence of interactions between fermions. Due to the 
interaction g 2 with the counterpropagating states, the re- 
flection amplitude will appear, being proportional to 32, 
thus the reflection probability is O (pf ) . This means that 
the transmission amplitude (and the probability) is of 
the form 1 — 0(g 2 ), i- e - the correction to the transmis- 
sion (and hence, to the conductance) vanishes to the first 
order in g 2 . 



2. Density of states renormalization: Diagram (e) 

The other fermionic self-energy that appears in dia- 
gram (e) in Fig. [2] arises from the 34 interaction. The left 
and right moving fermions interact only among them- 
selves via 34 which ensures that the resulting self-energy 
is translationally invariant, 



£ (e W) = - 9 g f' x>) ,v , 

llTl(X — X ) 



(3.22) 



as g4(x, x') = g±(x — x 1 ). This self-energy causes a per- 
turbative correction to the velocity of the scattering par- 
ticles. This is seen from the analysis of the Schrodingcr 
equation for the scattering states, 

-W»^.(x) + /^SW (a!> ^( a O=e0 e (a:). 

(3.23) 

The Schrodinger equation can be solved perturbatively 
with the help of the Wentzel-Kramcrs-Brillouin ansatz 



(3.24) 



ip e (x) cx exp [ia / dy ip{y) 



where in zeroth order tp^ — e/v. The first order correc- 
tion to the eikonal phase is determined by 
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In particular, the energy derivative of the eikonal phase 
up to first order is given by 

d l = '-^ (3.26) 

where 54 = J dx' ' g±[x — x')e % ^ x ~ x )l v is the Fourier 
transform of the g\ interaction at small wavevector e/v. 
This is equivalent to the renormalization of the velocity 
v — > v + 6v, where Sv/v = g±v. 

The velocity renormalization enters the normalization 
of the scattering states (|3.1a[> and (|3.1b[) and thus con- 
tributes to the conductivity by changing the density of 
states, v — > v{l — vg^). The density of states on the other 
hand enters the polarization matrix (|3.17[) which finally 
leads to the change in the conductance 

6G {e) /G LE = -2g A v (3.27) 

attributed to diagram (e). 



3. Random-phase approximation contribution: Diagram (f) 

The remaining RPA diagram (f ) in Fig. [5] can be re- 
duced to the spatial convolutions of the two components 
of the polarization matrix (|3 . 1 Of) . The resulting correc- 
tion to the current autocorrelator (|2.20|) reads 



G I 

K L G 









r < 






r = u 






r > 














— ► 



FIG. 3: Temperature dependence of two-terminal conduc- 
tance close to the Mott transition for different values of the 
tuning parameter r. On the metallic side, r < 0, at T — 
the conductance is determined by the Luttinger parameter of 
the leads, Kl\ Go = 2e 2 /h is the conductance quantum. In 
the insulating phase, r > 0, at lowest temperature the con- 
ductance is thermally activated with a universal exponential 
prefactor involving the fractional charge, eci = e/y2, of the 
critical degrees of freedom, Gci = 2eci//l- In the quantum 
critical regime, \r\ < T, of the phase diagram, see Fig. [TJ 
the conductance is determined by its critical value, G cr , see 
Eq. (|1.5j) . that depends both on the Luttinger parameter, Kl, 
and on the fractional charge, eci- 



4- Full perturbative correction 



5K, U) {x,x'- ) n n )= (3.28) 

- 92 e 2 cl v 2 J d Xl u^. :I , ir;! - n:', ;: <,. 11.: f :! > ; 

As we are only interested in the behavior at small fre- 
quency f2 n , we can use the simplified expressions (|3.17[) 
for the polarization operators. The integral over x\ can 
then be performed with the help of the convolution (no 
summation over index i implied) 

J dxiUl[ Xi . nn U^ i X ,. nn 

= ^ ; n„(A-^N-^l) (3-29) 

where (3 t was defined in Eq. (pTTTj) and EP-? in Eq. $3Q . 
After the evaluation of the x\ integral, we obtain for the 
RPA contribution 

n n ) = -2g 2 ve 2 Gl v 2 [-2v5{x - x') (3.30) 
— H-x,x';n n + n K) o;n n ^n 0)0; n n ^n 0) ,j,/.n,J 

The presence of the product of bosonic T-matrices, T, 
see Eq. (|3.18p . leads to the correction to conductance 
that contains an additional Fermi function, 

5GW/Gle - vg2 [2 - f(r/T)} . (3.31) 



Collecting all the perturbative contributions to con- 
ductance, (|3.21[) . (|3.27p and (|3.3ip . we obtain the follow- 
ing for the first order correction in the parameter Kl — h , 

6G/G = -^f(r/T) = (K L - |) f 2 (r/T), (3.32) 

where we used the identification vgi = 1 — 2K L valid 
to the first order in 52, cf. Eq. (|2.9p . As expected, the 
first order correction to conductance is independent of 
54. Note that the correction starts with the square of 
the Fermi function whose origin can be traced back to 
the RPA diagram (f) in Fig. [21 Adding to this the ze- 
roth order contribution (|3.20p we finally obtain the result 
announced in Eq. (|1.2p . 

IV. SUMMARY AND DISCUSSION 

In the present work we considered the finite- 
temperature dc conductance close to the commensurate- 
incommensurate Mott transition in a one-dimensional 
wire as a function of the tuning parameter, r, of the 
transition. For a homogeneous system in the absence of 
leads the result is universal and is given by Eq. (|2.2ip . 
It is determined by the critical degrees of freedom of the 
transition — the fermionic charge carriers that have a 
fractional charge eci = e/v2, — and the temperature 
dependence is simply a Fermi function. 
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Furthermore, we addressed the question on how the 
conductance is modified when the leads are attached to 
the Mott wire. The result is summarized in Fig. [3l In 
the incommensurate phase, r < 0, when the system is a 
metal, the dc conductance at zero temperature does not 
contain any information about the universal signatures 
of the Mott transition because it is solely determined by 
the compressibility of the electron liquid in the leads, 
i.e., by the Luttinger parameter Kl in agreement with 
Refs. and 0- However, in the insulating phase, r > 0, 
we found evidence that the critical degress of freedom 
of the quantum phase transition become manifest in the 
dc transport. In particular, in the lowest order in the 
fugacity e~ r / T , the conductance is independent of the 
properties of the leads, Kl, while the exponential pref- 
actor is determined by the universal characteristic of the 
Mott transition as it is given by an effective conductance 
quantum of the fractional charge carriers, Gci = 2e^j//i. 
The conductance in this regime originates from thermal 
activation of solitonic excitations that travel ballistically 
across the Mott-gapped region of the wire. Finally, in the 
quantum critical regime, |r| < T, of the phase diagram 
in Fig. Q] the conductance is dominated by its critical 
value (|1.5p that depends on the Luttinger parameter of 
the leads, Kl, as well as on the fractional charge eci- 

Our main finding is that the dc transport shows signa- 
tures of charge fractionalization (| 1 .4|) in the Mott wire. 
The underlying reason is that the electrons, irrespective 
of their interactions in the leads, must "transform" into 
Mott solitons with charge (|1 .4[) in order to overcome the 
Mott barrier A. The fractional charge can in principle 
be detected with the help of the finite temperature be- 
havior of the critical two-terminal conductance of a Mott 
wire. We contrast this fractionalization with that in a 
Luttinger-liquid wire, where only the forward scattering 
is included^ there, the phenomenon of fractionalization 
is not manifest in the dc conductance^^ 

The regime discussed above is restricted to the imme- 
diate vicinity of the quantum critical point of Fig. [TJ i.e. 
to the temperatures and to the tuning parameter much 
smaller than the Mott gap: \r\, T <C A. In particular, we 
neglected corrections arising from the thermal activation 
of particles in the other Mott band which are exponen- 
tially small in e~( 2A ~ r ^ T . Furthermore, we disregarded 
effects of tunneling of solitons through the Mott barrier, 
which gives rise to a finite conductance in the insulat- 
ing regime of the form SG ~ e~ s , where S cx rL/v with 
the length of the barrier L and the velocity of solitons v. 
The approximation of neglecting this tunneling correc- 
tion breaks down at sufficiently low temperatures when 
the length L becomes comparable to the thermal wave- 
length £_ T = v/T. 

We arrived at our conclusions by performing a con- 
trolled perturbative expansion around the special value 
Kl = 1/2 of the Luttinger parameter in the leads. At 
this specific Luther-Emery value, the solitonic excitations 
are free, and the problem of the Mott wire with leads re- 
duces to a ID scattering problem as descibed in Section 



IIII Al For Kl different from the Luther-Emery value, 
the interaction between the scattering states in the leads 
yields the result (fL2"|) . 

Our conclusions concerning the temperature depen- 
dence of the critical conductance differ from the results 
of the work of Mori et al.X There, a similar setup was 
considered, and the current autocorrelation function was 
also calculated in terms of the critical degrees of freedom 
of the Mott transition. However, the current correlator 
was only considered for a homogeneous system when the 
Mott gap extends across the full wire length resulting 
in a translationally invariant expression. The result was 
then used to match the current autocorrelator within the 
Mott wire with the one of the Fermi liquid leads in order 
to derive the conductance. We believe, however, that the 
use of the current autocorrelator for the homogeneous 
Mott wire is inconsistent with the matching conditions 
applied in Ref. 0, which explicitly break translationary 
invariance. Instead of using the homogeneous current au- 
tocorrelation function one may consider the result (|3.17|) 
for the case of a gapped region inside the wire as a start- 
ing point for a matching procedure. Preliminary consid- 
erations along these lines have led us to the conjecture 
described below. 

Let us speculate about the functional form of dc con- 
ductance beyond the leading Kl — \ correction. We con- 
jecture that, in the limit \r\,T <C A and T ^> v/L, the 
corresponding resistance is a sum of the two contribu- 
tions 



R — Rl + R 



A I 



with 



(4.1) 



h l-/(r/T) 



ci 



firm 



where eci = e/V2 and / is the Fermi function (|1.3| . The 
resistance Rl is attributed to the leads and can be inter- 
preted as a contact resistance. In particular, only Rl de- 
pends on the Luttinger parameter within the leads, Kl- 
The contribution Rm has the form of a four-terminal 
resistance with an effective transmission given by the 
Fermi function, f(r/T). It originates from the critical 
Mott part of the wire and involves the charge quantum 
of the critical degrees of freedom, eci- The expression 
(|4.ip reproduces the known limits. In the incommensu- 
rate phase, r < 0, at T = the resistance reduces to 
the contact resistance R — Rl- In the insulating state, 
r > 0, at lowest temperature the resistance is dominated 
by the contribution of the Mott region Rm and becomes 
independent of the lead properties. Expanding (|4.ip in 
Kl — \ one also reproduces our hard-to-guess perturba- 
tive result (fl~2ll . 

The Mott region in a quantum wire may be realized 
by employing an external periodic potential created by 
gating, 29 or an acoustic field.— 

Our approach applies also to the problem of drag be- 
tween two quantum wires. At equal electron densities, 
the inter-wire interaction locks the charges in both wires 
such that the drag resistance at T = is infinite in the 
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commensurate state. This system bears resemblance to 
the Mott insulator: if opposite biases are applied to the 
two wires, there is no current at T = . 31 > 32 Density im- 
balance drives the system into an incommensurate phase, 
for which one recovers the ballistic conductance for each 
of the wires; drag is absent at T = 0. The problem of a 
finite-T drag resistance may be mapped onto that of the 
critical conductance in a Mott wire. Technically, such a 
mapping is established by introducing the even and odd 
modes, p± — (pi ± p2)/y2, where pi^ix) are the elec- 
tron densities in the two wires. The dynamics of the odd 
mode is described by the sine-Gordon model (|2.ip , while 
the even mode is a conventional Luttinger liquid. The 
drag resistance is determined by the conductance in the 
odd sector. The present analysis of the critical behavior 
of the transport at the commensurate-incommensurate 
transition is fully applicable for the evaluation of the drag 
resistance. 
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APPENDIX A: POLARIZATION OPERATOR IN 
THE COORDINATE-FREQUENCY 
REPRESENTATION 

We present a detailed calculation of the polariza- 
tion operators (13.10[) . After performing the sum over 
fcrmionic Matsubara frequencies we obtain 



11 x,a;';n ? , 



dede' .^ (£) tr {t'A^,^ A a ,, a;e ,} 

iVL n - (e - e') 

(Al) 

with the spectral function A given in (|3.9p . In this sec- 
tion, we use a temperature dependent definition of the 
Fermi function, f(uj) — \j(e u l T + 1). As a first step we 
evaluate the trace. We will use the abbreviations (13. 8[) . 
Consider first a situation with s ^ s', i.e. with x and x' 
on opposite sides of the barrier. 

s ^ s' : tr {t 1 } (A2) 

= v 2 Q t - r Q c ^ r A l ^,{x -x' + sl t ) . 



Here 



(A3) 



where the 8 % were defined in (|3.11|) ; and we made use 
of the approximations (|3.6|) for the phase shifts with 
l t = (£+ +£-)/2. The case s = s' is done explicitly 



[in Eq. (|A4[) set v = 1 for brevity]: 
s = s : tf A. x x / - € t^ A. x f x - e '^ = 



(A4) 



v 2 tr 



a ie(x-x') 



s Q r _^ e ie(x+x')+is4> s (e) 



sf3 l Q r _ e e 



ie{x+x')—is(j) s {e) pi ^—ie(x-x') 

e —it'(x — x') S Q /e ie' (x+x')+istp s (e') 

sft Q r -e' e~ ie '( a: + a: ') -iS( ^( e ') fp g*^ 

= v 2 (A l l t ,{x - x 1 ) + B 3 Q r ^e r ^Al 3 el {x + x 1 + si.fj 

where we used the approximation (J3TBJ) for the reflection 
phase shifts. 

When substituting the above expressions for the trace 
into Eq. (|A1[) . it is convenient to define 



FJx) 



d£de m n - (e - e') ■ [ ' 



Note that in both cases, after evaluating the trace, the 
conditions e > r and e' > r, or e < r and e' < r, enter si- 
multaneously. Hence the polarization operator (|A1[) can 
be expressed via (|A5|) as 



Ll x,x'-,n n 
1 - ss' 
' 2 

1 + 6 



(A6) 



F r (x -x' + sl t ) + B l B 2 F r (x' -x- sl t )] 
\ 1L x,x'-n„ 



+ B' J [F_ r (x + x' + sl s ) + 8 l B3F- r (-x' -x- sl s )]\ . 

The last term is obtained from Eq. (|A5[) via the inter- 
change of the dummy variables e 



— e' and e' — * — e, 
and using /(— e) = 1 — /(e). We also define the polariza- 
tion operator in the translation-invariant case 



n 



(0)ij 



[F r (x - x') + l?FF r (x> - xj\ ^ . (A7) 



We now proceed to evaluating (|A5|) . Introducing the 
variable e = e — e', write 



F r {x) = v 2 



de 



^tex 



de 



/v de 



ifl n — £ 

■o „isx/v d£ 



/(e)-/(e-g) c i £a /^ 
iCln - e 

'de [/(e) -/(e- £ 



(A8) 



/■oo ie X /vj p roc 

/ V^7/ *[/(e)-/(e-e)]. 

JO W'n E Jr+£ 



' r+e 

The integration with respect to e and £ is over a domain 
that can be split into an infinite rectangle, e > r and 
e < 0, and an infinite triangle, e > r and < e < 
e — r. After changing the order of integration in each of 
the two sub-domains separately, such that the integral 
over e is performed first, one obtains the last equation in 
(IA8[) . Shifting variables in the energy integrations and 
introducing the function 



<p(e) 



r+e 



f(r + e) 

def(e)=s + Tln^-^, (A9) 
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we get 

F r {x) = V 2 



de ip(e) 



j—iex/v 



e — iQ n e + iQ, n 



(A10) 



Consider first the static part, F* tat (x) = F r (x)\ in =0 ; it 
can be represented as 



(x) = vf{r)8 B {x) 



(All) 



where Sb(x) is a function with a peak at x = that 
carries unit weight, J(1xSb(x) = 1. In the limit r —> 
— oo, the function i5b reduces to a delta function; in the 
classical limit, < r/T -C 1, the peak has a width of 
order £t = u/T. 

The dynamic part of the F r function, F^ yn = F r — 

^7>stat 



F^ n (x) = u 2 ifl n I de 
lo 



<P(e) 



a iex/v 



-iex/v 



e + iQr. 



(A12) 



is dominated by the poles, e = ±iSl„. The leading be- 
havior 



F^ n (x) ~ -7W 2 n n /(r) [signs + signf]„] e 



— |On3; I / v 



(A13) 



originates from the odd part <p(e) 



<P- 



ip(—e)] — ef[r) + C(e 3 ) by extending the integral over 
to the real axis and evaluating the Taylor-expanded 
<P-(e). The even part, <p + = |[y(e) + (f(—e)], and the 
higher-order terms of the odd part, yield the corrections 
0(f'(r)fl 2 l log that are subleading in the dc limit, 

uj/T -C 1, of small real frequencies, ifl n — * u> + iO after 
the analytic continuation. 



We now have all the ingredients to construct the polar- 
ization operator (|A6[) . First, the absolute values of the 
arguments of the F r functions in (|A6|) are all larger than 
the scattering lengths £± that characterize the length of 
the Mott region of the wire. In the limit of a long Mott 
wire, £±T/v > lwe can neglect the tails of the static part 
(|A11|) of the F r functions in (|A6j) . The static part then 
only enters the polarization via the contribution n^ ^-' 
which was defined in (|A7[) . The F r functions in (|A6[) thus 
only contribute with their dynamic part (IA13j) . More- 
over, for the dynamic part the scattering lengths £± in 
the arguments of the F r functions in (|A6|) can be effec- 
tively set to zero as they give only subleading corrections 
in the dc limit oj£±/v <C 1. After these simplifications we 
finally obtain the expression (|3.12[) for the polarization 
operator. 
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